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a topic requires specialized terminology. 


Problems designed to broaden and fix the concepts are supplied 
in the first seven chapters—algebra through calculus. 


Puzzlements are omitted in the interest of brief presentations of 
basic principles. 
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applications. 
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A cure for those who have the mistaken notion that they “just 
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for additional exercises). 


The Tree of Mathematics, containing 420 pages, 
with 85 cuts and pleasing format sells for the low price of 
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ABSTRACT STRUCTURE OF INEQUALITIES 


A. B. Soble 


INTRODUCTION 


Let f(z), g(x), and h(x) be any three single-valued rea] functions de- 
fined over the interval (z ” 2.) 

The reader may readily verify that A passes through all intersections 
(if any) of f and g in the interval, and otherwise lies between them, if and 
only if A is a member of the family of functions: 


Of + (1-49, 


where 0< O(x)< 1, (2, <a<z,). 
A similar result obtains if f, g, A, and 6 are functions of the two real 
variables z and y over a given region. 
In the remainder of this section only, let x+y and xy represent the 
l.u.b. and g.l.b., respectively, of the two functions 


¢ 


More abstract concepts of betweenness are as fol 


Let z, y, 2 be distinct functions, with y distinct from 2+2 and from a 
Then any of the following relations is a necessary and sufficient con- 


dition that y lie between z and 2, except at their intersections. 
yYln+ 2) =Y=Y + x2 
ry + YZ = ¥: (x4 yy + 2) 
G+ 22Y 2 %2 


a+yz2>y>xly+ 2). (4] 


We note that the conditions are self-dual, since interchange of “addi- 
tion” and “multiplication”, and reversal of direction, leave the conditions 
unchanged. 

_Betweenness is an order relation of three elements. In this paper we 
shall consider the order relation of two elements, called inequality or 


majorization. 
ONE DIMENSION 


Let f(z) be a single-valued real function defined for all real z. 


Let z, < By <u S By < By. 


Suppose it is desired to find the family of real functions g(x) such that 
179 
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in the closed intervals (z,, 2,), (%,, %,), and (z,, 2), g(x) > f(x); while 


outside these intervals g(x) < f(x). 
It may be readily verified that the solution is 


g(x) = f(x) - G(z)Il (x-2,), 
=1 


where G(x) is any real function which is positive outside the intervals, 
and positive or zero inside. An example is G(z)= 1. 
Extension obtains for any number of non-overlapping intervals. 


TWO DIMENSIONS 


Let f(z, y), F (2, y), (é=1, 2,3), be single-valued real functions defined 
over the real plane. 

Let F (x,y) = 0, (é=1, 2,3), denote the boundaries of three closed re- 
gions which do not overlap. 

Since the boundaries F’; = 0 are the same as -F’; = 0, we may assume 
that the signs are chosen so that the interiors of the regions are given by 
F,<0. 

Suppose it is desired to find the family of real functions g(z,y) such 
that in the closed regions F(z, y)< 0, g(x, y) > f(x, y); while outside these 
regions g(z,y) < f(x, y). 

It may be readily verified that the solution is 


3 
g(x,y) = f(z, y)- G(z,y) MF, 


where G(z, y) is any real function which is positive outside the regions, 
and positive or zero inside. An example is G(z, y) = 1. 

Again, the result may be extended to any number of non-overlapping 
regions. 

There remains the question of obtaining the equation of a region when 
its boundary is given by a graph, map or numerical table. 

If the region is star-shaped, that is, there exists an interior point 
(z.,¥Y_) such that every ray emanating from (z,,y,) intersects the bound- 
ary in one and onlyone finite point, the equation may be found as follows: 

Let the angle of the ray be 6, and let the distance from (z,,y,) to the 
boundary be p. 

Since the numerical value of p is known for -7 < @ < 7, it may be ex- 
panded in a Fourier series. 

Hence p = p(sin@, cos§). 

Therefore, 
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Y-Yo TZ oy 
V(e-zx 7 +(y-y,)? = , 

0 0 y) 3 a 2 
(r—-z,) + (y-y,) V (2-2) + (y-¥9) 











is the equation of the boundary. 


TRANSFORMATIONS 


Let F(t) represent the class of real single valued functions defined for 
all ¢ in the half-open interval ¢, < ¢ < ¢,, and lying between the constant 
bounds « and # as in Figure 1. 


nme fp 





F(t) 


Pe 











Figure 1 
7 
The transformation 6 = ———[(¢-t,)+(¢-¢,)] and p = F-« carries F(t) 
bom 
1 
into a star-shaped curve about the center of, and lying within, the circle 
of radius B—«, as in Figure 2. 





“ed 





Figure 2 


If G(t) > H(t) in the rectangle of Figure 1, then the transform 


surrounds or contains the transform of H in the circle of Figure 2, 


conversely. 
Let g (F) be any monotonically increasing function of F, « < F 





182 MATHEMATICS MAGAZINE (March-April 


Let y and 6 be constants, and S(¢) a function, such that y < %(t) < 4, 
when ¢, <¢< ¢,. 
If y + @(«) > « and 5 + (f) < B, then S(t) + @(F) is a majorant-pre - 


serving transformation of the class of curves F(t) into a sub-class. 
In this paper we shall consider algebraic properties of the transfor- 
mation 


F-« «+ BF 
or -_——+ 


* B-F 1+F 


7 


which converts the bounds («, 8) into (0, ). 
Since 
OF B-« 
ae eee > 
OF (B-F)? 
e. increases strictly with F. 
Hence the transformation preserves majorization. 
We note that the transformation gives a 1 to 1 correspondence between 
F and F 


We shall also consider algebraic properties of the transformation 


F- —«F 
»* B or F =<. 
PF —« 1-?, 


0, 


F 


which converts the bounds («, 8) into (—s0, 0). 
Since 
ae salle tey > 
OF (F-«)? 
increases strictly with F. 
Hence the transformation preserves majorization. 
We note that the transformation gives a 1 to 1 correspondence between 


F and Fo: 
OPERATIONS 


Let G(t) belong to the same class as F(t). 
We define the commutative binary operation F ® G as 


«+B(F +G) 
140. #G.) 
” 7 


0 


, 


F 


0 





Then (FOG), = r. + G 
Since F andG >O,then(F @G) >F_. andG.. 
+ *x*~ = 7 * 


Because of the 1 to 1 correspondence and the monotonicity, it follows 
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that F © G >F andG. 
It remains to show that F © @ is also of the class of F. 


Now 
Paw G -—« 


v J 
p-F 6-¢ 
is defined and single valued whenever F and G are. 


That F © G is also, follows then from the definition of F 6 G. 
We must still show that«<F OG < B. 


Since F 6 G > F and G, it majorizes «, 

Moreover (F © G)«e = Fe+ Ge <o. 

Hence, by the 1 to 1 correspondence and the monotonicity, F @ G < B. 
Thus F © G is in the original class. 


(Fe + Ge) = 


We define the commutative binary operation F 0G as 


B-«(F,+G,) 
1-(F,+G,) — 
Then (F OG), = F, + Go. 
Since F and G, < 0, then (F OG), < F, andG). 
Because of the 1 to 1 correspondence and the monotonicity, it follows 
that F OG < F and G. 
It remains to show that F © G is in the original class. This we leave 


to the reader. 
ALGEBRA 





We have already seen that the commutative binary operators 9 and © 
are operators of majorization and minorization, respectively, namely 


'FO@G>FandG; FOGS<F andG. 


It may be readily verified that « and f are the identity elements of ® 
and © respectively, namely 


«@F=F and BOF=F. 


There is no inverse. 
Since 


(FOG) OH). = Fa+Ge+Hs = [F O(G6H))s , 


it follows that the operator © is associative. 
Similarly the operator © is also associative. 
The two operators are not mutually distributive. 
We define the complement F’ of F by «+B-F. 
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The reader may readily verify that: 

(F’)’ = F, «= B, and B’ = « 

(F’), = -F, and (@°), = -Fe. 
"e(F’)y = FCF’), - and FFs = -1 
(F @G)’=F’OG’ and (F OG)=F’°OG. 
The algebra which we have defined we may call an additive algebra, 
because 
(F ® Ghe = Fe + Ge 


(F 0G), =F, +G 


0 0° 


A multiplicative algebra may be developed by defining (F © G)* = F*G* 


and (F © G)° = F°G°, under the transformations 


(«, B) + (1, 0) 


(«, 8B) + (0,1). 


This we leave to the reader. 
Still another algebra of majorants is that given by lattices [1], in terms 
of l.u.b. and g.1.b. 
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THE COEFFICIENTS OF £2S42 AND A NOTE ON 
CARLITZ’S COEFFICIENTS OF Sithz 


sin x 


J. M. Gandhi 


INTRODUCTION. The present paper is divided into two parts. The 


f cosh z 
cos 


=< and of sinh x 


sine will be discussed 


properties of the coefficients o 
separately in two parts. 
PART 1 
THE COEFFICIENTS OF cosh z 


xr 


a7 2n 


cosh z _ s nd 
cos = 2n (Qn)! 
R= 
“es x? na 2 2n vi 22 
= 42, (~1)" S Ss, —— |. 
(2n)!] — (2n)! 2n (On)! 


n=0 n=0 


Whence 


Simplifying and equating the coefficients of 2?” 


- 9 
an 
) -1)Y : 
ie a rs 


y=0 
From [3] we get : 


we get 


2n Qn [on 
anf 2” \on- Son—g t ++ (-1) S,=1. 


Solving for S,,, and using Cramer’s rule we get 


185 
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] 0 
~(Z) (3) 


()) ~(4) 





0° (<n (% ... (92%) 





TABLE 


S,=2 S,=12 S,=152 


S,=3,472 5S, = 126,752 


Now we shall prove that 
(A) All S’s are even positive integers. 
(B) The last digit of S,, is 2 for n > 1, whence S,, - S,,_,, =0 
(mod 10) for k < n. 
(C) S,,, = 0 (mod 4n-1), if 4n-1 is prime. 
Theorem (C) can be easily proved from formula [3]. Proof of (A) : re- 
write formula [3] as : 


Bi, a ag be on gt 20 =U, +1, [5] 


If S,, Sy, --- Son_a are integers we see that S,,, is also an integer, but 


from the table we see that S,, S,, ... S,,_., are integers for some value of 
n, hence they are integers for n+1, whence the theorem follows by induc- 
tion. 
Also if S,, S,, ... S,,_, are even and since (-1)°(2%)5,41 is 2 or zero, 
S,,, are even and the theorem that all S’s are even follows by induction. 
We shall prove [See page 5, formula 18] that S,,, = [EF +1]?", where E’s 
are Euler’s numbers and all are positive whence S,, are all positive. 


Proof of (B) : rewrite formula [3] as 
Ban = gS 55 oboe Sanna too + -1)S, +1 
an ~ ‘2n—2’" an—27 ‘2n—4’" an—4 "°°" “Voit 


Let us first consider the case when n is odd. 
Then 


2 2 2 
Son = (on-2)Sana~ (ond Sanat ++ (Z)Sy + 2 


Assume that 
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S,, S,, 8, --- Sa, = 2 (Mod 10) 


2n 


Then from [6] we find that 


Say = 2+ 252%) - (924) +... (2) - (29)] [Mod 10] 


Since the value of the bracketed term is zero we get 


San=2 (Mod 10) [8] 
From the table we see that if [7] is true for some value of n, it is true 
for n+1, whence the result follows by induction, provided that n is odd. 
Now consider the case when n is even. This case is an interesting one 
as here are used the properties of the imaginary quantity i, i.e. V-1. 
Let n = 2k, formula [3] becomes 


4k 4k k 
San = (4p )S cpa apra Sapna tee + QS, [9] 


Assume that : 


S., S,, --- Sg, 9 = 2 (Mod 10) [10] 


ee 


Then from [9] we get : 
Say = Ug po) - ghk,) +...+ (41 [Mod 101 1] 


By using Binomial theorem the following identity can be easily proved. 


(142) 4"+(1-1) *” 
4k 4k 4k), (4k) _ 
(gpa) — (gpg) ++» C4) + C9) = 2- 9 [12] 


a 





Since (1+¢)? = 2¢ and (1-1)? = -2i, formula [12] becomes : 
(fF) - (ht) +... OR) + (HH) = 2-(-1)*427* 
Substituting [13] in [11] we get 
S,, = 212-(-1)*2?*] [Mod 10] 
Consider first that k is even. Let k = 2m. Then 


S 4—2(16)” (Mod 10) 


8m 


= 4~2.6 = -8 = +2 (Mod 10) 


Now consider that & is odd. Let & = 2m+1. Then 
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4+8(16)” (Mod 10) 


4+8-6 = 2 (Mod 10) [16] 


Now from the table we see that [10] is true for some value of &, from [15] 
and [16] it follows that it is true for k+1, whence the theorem follows by 
induction. 

From the table we make the following conjectures: 


(A) n(2n-1)8,,,_ 5 > So, 


(B) S,,,= 0 (Mod 2”). 
Conjecture (B) can be proved if we can prove that k,, are positive odd 
are given by 
oo 
cosh—— 
on KS 
cos 


V2 n=0 


integers, where k, 


Now we shall try to construct a formula relating S’s with the famous Euler 


numbers. 
Let 


oo 


2np 

l ; ; 7 Fad sii 

cos z (2n)! 
n=0 


With this definition of Euler’s numbers, all of them are positive. 


Now from [1] and [17] we get : 
2 


~~ et. 2°” 
> ; > Gay 2] - > - 
(2n)N] 2, (2n); & (2n)! 


n=( n=0 


fo] 


92 
S i" 


Equating the coefficients of x*” and after some manipulation we get 
S,, = (E +1)". [18] 
Where after expansion, the index is to be replaced by the subscript. 


Equation [18] gives a method to calculate S’s if E’s are known. Inci- 


dently [1] also implies that 


=e 2n 
cos z 
oo. - S (-1)"S.,, 
cosh z (2n)! 
n=0 
so that from [1] and [19] we get : 


(-1)” (5%) s S,,, = 0. 


an—2y2y ~ 
y= 





COEFFICIENTS OF S¢882 


PART Il. 


A NOTE ON CARLITZ’S COEFFICIENTS OF 


Let 
° oo an 
sinh z ‘ x 
* (Qn)! 
n=0 
Whence we get : 
= ns go 2M 1 


yg ent} si pen 
SS nS Se, 
4 (n+)! (2n+1)! 7 Coa! 
n= 


N= 0 n= 0 


Equating the coefficients of r?"*! and after some manipulation we get 


nr 


> (DY GOD Ban ay = 1 [23] 

y=0 5 
As done in part 1, from formula [23], we can also obtain a determinant 
for B’s. The determinant for 6’s had already appeared in MATHEMATICS 
MAGAZINE as a problem (See reference number 2), 
Formula [23] gives an easy method to calculate £’s. Professor Carlitz 

has calculated £’s from the formula 
m 
Imaal , 
Bin =e ZO (“ar")D,, [24] 


This formula for calculating 8’s is tedious, as for calculating B’s we 


must know D’s, while formula [3] directly gives the values of B’s. 


TABLE 


2 _ 3104 79,808 
Prop Pia 


It is evident that Professor Carlitz’s table for 8, and A, is wrong or 
misprinted. 
We shall prove that : 
(A) All B’s are positive. 
We shall prove that : [See formula 28, page 191] 
(B Pe ace = at Por 


We have already proved that all S’s are positive and if B,, B,, ... By,» 
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are positive then B,,, are also positive. Since B,, B,, ... By,» are pos- 
itive for some value of m, it is true for m+1, whereby the theorem is es- 
tablished by induction. 
(B) Theorem : If 4n+1 is prime, em will never contain a factor 4n+1 
in its denominator. (A special case of Carlitz’s theorem, but here is a 
simple proof for it.) 
For from [23] we get 
2n-1 
An+] 
(-1)” (Oy4D) Bin-ay = 0, 
y=0 


2n-1 
(4n+1) By, ao > (tp) Ban-ay’ [25] 
y=1 
Since when 4n+1 is prime, the right hand side of [25] is exactly divis- 


ible by 4n+1, and since B,,,, when 2m < 4n, will not contain any factor > 
(2m+1) in its denominator, hence the result. 
From the table we make the following conjecture. 
(1) Numerator of 8, = 0 (Mod 2”). 


Now rewrite equation [21] in the form 





(sinh. cosh> 3) S tuck 


(sing 5° cos 3) m=-0 


2m (2m)! 


Whence we get 


~ am 2m . 2m 
> Pom ser) a il > io 7m = Percy 
m=0 m=0 0 


Equating the coefficients of 2?” and after some manipulation we get 





n 
(B+1)?” = > aie (-1)”. [26] 
y=0 
Now we shall construct a formula relating S’s with f’s. 
Formula [21] can be written as : 


(sinn 5. cosh> 5) >, B 
2m 


(sin5 2° *COos 9 =) m 





oer 
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Using [1] and [21] we get from [27] 


sie 22m P 22m 22m 
> Pam 92™ (9m)! > am S7™Im)!| Pam (2m)! ° 
m=0 m=0 m=0 ee: 


Equating the coefficients of 27”, and after some work we get 


(B+G)?” = 29" 8 [28] 


2m°* 
Where as usual after expanding, the index is to be replaced by the sub- 
script. 

It is interesting to note that similar results hold good for various hy- 
perbolic, trigonometrical and inverse functions namely coefficients of 


tanh z coth z sechz 
tanz’ cotz’ 


nominators in several cases possess the same properties as that of 


~ = , shot 
etc. and Simh”_2@ cosh” 2 ote, and that the de- 


, 


sec z sin 2 cos 2 








Carlitz’s coefficients. 
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ASME TO TRANSLATE , RUSSIAN JOURNAL 


A leading Russiantechnical journal will °°? be translated into English 
of Mechanical Engineers, it 


on a regular basis by The American Society ~ 
was announced recently. The effurt marks the fi, "St time that the engineer- 
ing society has undertaken such a project. : 

Under a $35,000 grant from the National Science b oundation ASME will 
publish the bi-monthly “Journal of Applied Mathematics and Mechanics.” 

The announcement noted that the Society had underta.ken the task of 
translation “in an attempt to correct the present situation !n which the 
Russians are familiar with the content of most, if not all, of our technical 
publications, while only a few of theirs are translated for use by the 
English-speaking world.” 

Arrangements have been made with Russian scientists, said an ASME 
spokesman, to secure proof shééis of the Russian journal in advance of 
final printing, to permit speedier translation. 

In the Russian language the magazine is known as “Prikladnaya Mate- 
matika i Mekhanika”, usually abbreviated as “PMM”. 

The magazine contains the latest theoretical and practical advances 
made by Russian scientists in mathematics, fluid dynamics and solid state 
physics. Copies will be sold, on a subscription basis, to any interested 
persons or groups at an annual rate of $35 for the six issues. ASME mem- 
bers are entitled to a 20% discount. Subscriptions may be ordered from the 
Order Department, The American Society of Mechanical Engineers, 29 
West 39th Street, New York 18, New York. 

In a statement accompanying the announcement, James N. Landis, 
ASME President said, “Translation of the Russian Publication ‘Journal of 
Applied Mathematics and Mechanics,’ one of the leading publications of 
its kind anywhere in the world, is expected to make important technical 
contributions to the English-speaking world. It will provide many English- 
speaking eagineers and scientists with access to the latest U.S.S.R.de- 
velopments in a truly fundamental field. 

“In addition to making information available, however, the very initia- 
tion of such a project is a valuable contribution in itself. It serves to un- 
derline the universal nature of science and engineering, and to emphasize 
that in these fields there are no natural boundaries between nations. It is 
my hope that this project, together with the benefits that are sure to follow, 
will help to break down some of the artificial barriers that now exist, and 
that from such relatively modest beginnings as these, better international 
understanding will grow and flourish.” 

Professor George Herrmann of Columbia University will serve as editor 
of the translated magazine. 





THE MIDPOINT METHOD OF NUMERICAL 
INTEGRATION 


Preston C. Hammer 


It is a noteworthy fact that while virtually all calculus texts and nu- 
merical analysis books discuss the trapezoidal method of numerical inte- 
gration they ignore one which is simpler to use and generally superior to 
the trapezoidal—namely the midpoint rectangular method. The midpoint 
method also may be used with the trapezoidal method to obtain bounds for 
an integral in some cases. The midpoint method is an application of the 
Newton-Cotes “open” formulas with one point, 

The trapezoidal and midpoint methods over one interval of length / 


given respectively in the following equations : 


an 
1) | fla)de = 5(f(. 
Zo 


oo 


I 
7 J f(x)dz / f(x, i * 
xr 4 


0 


In these equations, the second derivative of the integrand is assumed to 
exist on [z,,z,]. It is clear that the error in the midpoint formula is about 


% the value of the error of the trapezoidal as z,+z, if /” is continuous on 
(z,, 2 ,].However, we now prove that 2) is superior to 1) always if f~ does 
not change signs on [z,, z,]. 


THEOREM: Let f be a function either concave* or convex, and continuous 


on [z,,2,]. Then 
1 h 1 
f(a) dx -5 (te) + f(a, ) >| (  f(a)de - hf(e,,) 
Zo ry vo 


and equality holds if and only if the graph of f(x) is a line or a broken 
line with break point at z = z,,. 





*A function is concave if any line segment joining two points on its graph 
lies below or on the graph. 
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‘ 
Proof: If {” exists, which we do not assume, and f is concave, then f”< 0 
and if f is convex /”30. The argument for either case being the same, we 
will assume f is concave and since adding a constant function to f does 
not change the error for either formula we may assume f >0Oon the interval 


[z,, 2,1]. 














7 


Figure 1 





Consider the graph in Figure 1. Since f is concave, the area bounded 
by y=0, z=2), r=2,, and y=ffz) is a convex set. Now the midpoint for- 
mula gives the area under any line through P,, over [z,,2z,] (more exactly, 
the integral of any such linear function), Choosing m, as a supporting or 


tangent line to the curve at P,, we have 
1 
A,+A, 2 Mte,)~[ f(z)dz >0 
Zo: 


where A, = area PoPyQ, and A, = area Py P,Q ,- Now, on the other hand, 
the trapezoidal formula gives an approximation which is less than the 


integral and 


| h 
J f(a)dz -(f(2,)4 f(a ,))> A, = area P,P Py 
x 
0 


But if we now take line m, through P,P, we readily see A,=A,+A, = area 


P.Q.P, = As;. Hence 








NUMERICAL INTEGRATION 


-, h ai 
f(a) dz -5((2,)+ 1@,) > Mt) ~ | Ka)dz. 


a) “ 


J 


Equality can happen only if f is linear or if the graph of f consists of 
two line segments joined at P,,. Hence, the midpoint formula is superior 
to the trapezoidal for either concave or convex functions and has an error 


opposite in sign. 


COROLLARY : If f is concave and continuous on [x , 2 | then 
nm 


n~| z, n h 
3) h f(x;, ra f(akde>h) f(x;)-—f(z,)+ f(z,) 
Y ey? | > fey -Zfleo)+ Me, 
i=0 ' i=0 
and 


™ 
fla )dzx 


x 
Zo 


n-1 ” 
4a) |r + 3 fleinyg)-[ "Me)de < IS le) fie, fle) -| 
i=] 79 0 : | 


where the equality signs in the latter case hold only if the graph of fis a 
line or a broken line with vertices at P;,,,. If f is convex the inequality 
4) holds and the inequality signs are reversed in 3). 

Simpson’s rule results if we take a weighted average of the midpoint 
and trapezoidal rules giving weight 2/3 to the former and 1/3 to the latter 
That the midpoint rule is easier to use than the trapezoidal in genera 
either for hand calculation or by machine, is obvious unless there 


reasons why f(t 544) is more difficult to compute. Calculating bot 


case f is either convex or concave, would provide an error bound for the 
integral. 
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A DETERMINANT FORMULA FOR HIGHER ORDER 
APPROXIMATION OF ROOTS 


J. M. Wolfe 


1. Introduction. Newton’s method of approximation of a root of a function 
replaces the original function y= f(z) by a linear function having the same 
value and the same value of the first derivative as the given function at a 
value of z near the root sought. In this article a convenient determinant for- 
mula is presented for the approximation when the agreementof the approx- 
imating function with the original function is extended toderivatives of 
higher orders at the point of contact. This formula is based on the method 
shown by Halley (1694) [1], Lagrange [6], Willers [11], and many others, 
for example, Stewart [9] and Snyder [8]. The method has been extended to 
n equations in n variables by Frank [4], for example. 

2. The approximating function. Let the approximating function be 


2 n—1 n 
(1) P= Ay + Ayr agy +e +G,_ yy +a,y 


where the values of the coefficients can be determined by taking z=2, 
and y and its derivatives equal to the corresponding values derived from 
the original function at the specifically selected point (z,,y,) near the 


root sought. 
The simultaneous equations determining the coefficients in the approx- 
imating function are 


n—1 4) ay 
[a + Ya, torr + tae t y a, = 


0 + Diya, +o + Diy" _ + D(y")a,, 


0+ D*ya,+--+ Dy Ya, + D%*y"a, 


0 +D™ Myla, +--+ +D™ My™ a, ,+D™ Uy")a, = 0 





LO + D™y)a, +--+ D™y™a,_,+ D™y"a, = 0 


where z=2,, y=y,, and D*(y/) represents the value of the i-th derivative 
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of ys with respect to x at the specifically selected point (z,,y,) near the 
root sought. 
8. The formula. When y = 0 in the approximating function (1), 2 = .. 


Thus it is necessary to solve the system of equations (2) only for a), 
whence 


2 eee n—1 n 


y y y y 
v*(y) D*(y?) D*(y"—") D*(y") 


Dp" 'y) Dp" iy?) p"™- iy" ‘s p*- iy”) 
D"(y) D'(y?) eee D'™(y"—') D™(y") 
D(y) D(y?) eee Diy") Dy”) 
D*(y) D*(y?) D*(y"—') D*(y") 





D™*'(y) D™ "(y?) oa *) Dp" My") 
D"(y) Dy?) see D'(y"—') D"(y") 








provided that the denominator is not zero. 

The following table demonstrates the efficiency of formula (3), taking 
n=4, It consists of three problemssin which Frame [2] used a formula by 
Halley [1] with a correction for agreement as far as the second derivative. 





Function y=2°-2 y = 82° -62+1 y = @-cose 





Root 1.25992105 . 766044 
(decimal approximation) 





at 








Error by the Newton- 
Raphson formula 

J ‘ 
r=2, ed .00007895 .000623 .000451 


, 


Y 1 
Error by a formula by 
Halley 





-.00000042 ~.000018 .000011 





Error by formula (3), 
taking n=4 -00000000 -000002 .000000 

















APPROXIMATION OF ROOT: 
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SYRACUSE PROFESSOR NAMED HEAD OF YESHIVA 
U. MATH INSTITUTE 


Dr. Abe Gelbart, a Syracuse University professor, who was inspired to 
pursue a life-long career in mathematics through a chance meeting with 
the late Professor Jekuthiel Ginsburg in New York’s Public Library at the 
age of 11, today succeeded his mentor as Director of Yeshiva University’s 
Institute of Mathematics. 

An expert in the fields of analysis and applied mathematics, and co- 
developer of the important theory of pseudoanalytic functions, Dr. Gelbart 
also will serve as editor of the University’s noted scholarly journal, 
Scripta Mathematica, which Dr. Ginsburg founded and edited for twenty- 
five years. 

Dr. Gelbart, who will hold the rank of professor of mathematics, fills a 
post left vacant since Dr. Ginsburg died last October at the age of 68. 

Possessor of a fervent interest in mathematics from his earliest years, 
Dr. Gelbart first met Dr. Ginsburg in 1922 when he was 11 years old. A 
resident of Paterson, New Jersey, he was periodically brought into New 
York City by his parents for visits to the New York Public Library, at 
Fifth Avenue and 42nd Street. There, while each was browsing through 
mathematics books, he was noticed by Dr. Ginsburg. 

Fascinated by the zeal of the scholarly youth, Dr. Ginsburg arranged 
with the boy’s parents for weekly visits to his office and home in New 
York. “There we would engage in informal talks about a great variety of 
subjects, actually anything I wanted to discuss, but especially about 
mathematics,” Dr. Gelbart recalls. “He kept me stimulated, always intro- 
ducing new and exciting facets of the subject. But most important he 
guided my reading, insisting that I read only those works which were sig- 
nificant. In that way, he greatly accelerated my learning.” 

During subsequent summers, Dr. Ginsburg invited the youth to his 
cottage in the Catskills. Their chance meeting grew into a warm friend- 
ship, lasting until Ginsburg’s death. 





KNOCKING A CONE INTO 
A COCKED HAT 


Daniel B. Lloyd 


A. few eyes, ears and horns attached tothe accompanying figure (Figure 
1) would readily transform it into a frightful creature indeed. But the author 
is content to leave it thus unadorned and merely let the reader imagine 
such unsavory possibilities. It is not considered good technique to dis- 
courage the reader too early in the article, much less to scare him away 
by the first picture. Just to dispel any vestige of naivité that “there aint 
no such animal”, it might be of passing interest to mention that a certain 
well known toy manufacturer finds it profitable to turn these out by the 
thousands, and much to the delight of our younger generation. 


Vo 


Figure 1 

Of more than passing interest, however, is the mathematics involved in 
designing these creatures. One has a right elliptical cone and cuts it by 
a series of planes inclined to the axis of the cone at such an angle @ as 
to produce circular sections. If a rotation of 180° is then made at each 
joint, an accumulated bending of the axis results to produce the distortion 
shown. Since each joint is circular, they turn freely on each other but are 
not readily detachable, as the rims are designed so as to fit neatly into 
one another. 

This property finds usefulness in the design of sheet metal flues such 
as are found in intake and offtake gas lines associated with smelting 
plants and large chemical factories. Here it is helpful to have the bend 
joints of circular cross-section although the flue itself may be elliptical. 

We desire to find the angle that the cutting planes make with the axis 
of the elliptical cone, in order to cut circular cross sections. We propose 
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to find a general expression for the angle @ in terms of « and B, where 2« 
is the least angle and 28 is the maximum angle at the vertex V of the right 
elliptical cone. Figure 2 shows front, side and top views, and also an 
auxiliary orthcgraphic projection of one of the oblique circular sections. 








Figure 2 


Referring to the side view, we have at once, 


Pi a-k 
sin /2 z sin /3 
£2 = 90°+«, (£3 = 90°—(«+6) 
(a-k) cos « 
| =———— 





a 
cos (« + 4) 


lr. a+k 


sin/l : sing 





{1 =90°-«, /¢ = 90°+(«~6) 


(a+k) cos « 
elo Sse 


cos («— @) 
If the oblique section shown is circular, r, =r, = 
(a*-k*) cos? « 


i cos (« + 9) cos («- 6) 


= 7" 





The locus of the right section, in its XY-plane is: 

2? y? 

or ee | (2) 
ae 38 


Let L (=-k,r) be one of its two intersections with the circular section. 
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Substituting it into (2): 


re » (a7-k*)b? 
32+ 527) and +f edn cad 


But a=/tan«, and b=hAtanB 


is (a?-k*) tan? B 


tan? « 


2 





a 


Equating (1) and (4), 


cos? « tan? B 





cos («+0) cos («—6) ~ tan? « 


cos? « tan? B 





Y(cos 2«+cos 20) tan?« 


2sin?« 


cos 2« + cos 26 = . 
tan~ 8 


2sin?« a 
cos 26 = 77 +2sin*«-1 
an“ B 


‘ sin? « 
cos* @ = 


aa | Eo 2 
+sin* « = sin“ «(1+cot* 8) 
tan? B 


6 Sin« B 
.0os0 = >« 
sinB 
Another solution, interesting by comparison, is the following: 
With V as the origin of cartesian coordinates, the cone’s equation is: 
2 2 2 
z 2 


a? . 6? 2 4 (5) 
Rotating axes through the desired angle 6, 
z= 2,cos@-z,sin@ 
Y=V, 
2=2,s8in@+2,cos0 
Then (5) becomes : 


cos? sin? 2 ysin?@ cos? 
(SS - ar pete 73 ea "Na8 + (ou) 22 ,=0 (6) 








For 2, any constant (40), we ask that (6) shall yield circular sections. 
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This requires: 
q : cos?@ sin?@ 1 
a? rs eee 5? 





Recalling (3) above, we have finally: 
Sin« 
cos@=— (8) 
sinB 

Evidently @ increases with 8 -«; or, the greater the eccentricity of the 
right section, the more exaggerated will be the resulting bend in the axis 
of the cone. For a right circular cone, 6=0, and no bending is possible. 
The results are general for « and 6 within the open intervals 0<«<90°, 
0<B<90°, 

The amount of bending of the axis of the cone is 26 for each cut, so 
that any desired bending may be obtained, — within limits. The size of the 
maximum vertex angle of the cone will determine how many cuts can be 
made without them overlapping. 

Evidently, if 6>90°-« the section cut would be a hyperbola. The limi- 
tation 6<90°-« is also apparent from formula (8), wherein if 0=90°-«, 
B =90°. 

More general cases of rotation,—other than 180°, lead to interesting 
possibilities which we shall leave to the reader’s investigation. We close 
with the observation that the cone, as a never-failing source of interest, 
has been a favorite lure of mathematicians through the ages. 
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NORMAL CURVE AREAS AND GEOMETRIC 
TRANSFORMATIONS 


David Gans 


It is customary to enploy the methods of the calculus in order to obtain 
the statistically important knowledge concerning the areas under the nor- 
mal curves, i.e., the curves of the family 


1 * 
I = e~% (2-2 )? | 
where z and o are parameters taking on all real values, excluding o=0. 
Viewed -broadly, this knowledge .can be regarded as consisting of two 
facts : 


(1) the tabulated values of the so-called normal curve areas, i.e., the 
areas under the particular normal curve (o=1, 7 =0) 


2 
et . 1 


] 
(II) met. 

(2) the fact that the area under an arbitrary normal curve (I) from 2, to 
z, is precisely equal to the area under the curve (II) from ¢, to ¢,, where 
t, = (@,-2%)/o and ¢, = (2, -%)/o. 

While there would appear to be no alternative to using the methods of 
calculus in order to obtain the knowledge in item (1), the same cannot be 
said for that in item (2). In fact, as will now be shown, this knowledge can 
also be obtained, and perhaps more easily, by using familiar properties of 
the simplest geometric transformations. 

Thinking of equation (I) as representing a specific curve corresponding 
to any fixed values of Z and o other than 0 and 1, respectively, let us- 
subject the plane successively to the following three transformations : 


(a) 2’°=2-%, y’=y; 
(b) ig r/o, y= y's 


oe 


(c) f= 2", a =oy"’. 


These transformations represent, respectively, a translation, a horizontal 


1For convience we use ¢ instead of z. 
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one-dimensional strain, and a vertical one-dimensional strain. The re- 
sultant, 7, of these transformations, namely, 


(III) t=(4-2%)/o, 2= oy, 


is a continuous (.ffine) transformation which, as is easily verified, sends 
the curve (I) into the curve (II). Also, T transforms the locus y=0 into the 
locus z=0, and z=a constant sinto ¢=the constant (k-%)/o, that is, the 
horizontal axis goes into itself and vertical lines go into vertical lines. 
Since the transformations (a), (d),(c), respectively, preserve area, multiply 
area by 1/o, and multiply area by o, the transformation 7 preserves area. 


y (or 2) 





It follows th «t the figure ABCD, where CD is any arc of (I), AC, BD are 
the ordinates of C, D, and AB is a segment of the horizontal axis, will be 
transformed by T into the figure ABCD’, where CD’ is the correspond- 
ing arc of (II), AC’, BD’ are the ordinates of C’, D’, and A’B’ is a seg- 
ment of the horizontal axis, and hence that the areas enclosed by the two 
figures are equal. If the abscissas of C, D are z,, x, and those of C’’, D’” 


are ¢,, ¢,, then, according to (III), ¢, = (2,-%)/o and ¢, = (x7, -%)/o. 
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COMBINATORIAL DERIVATIONS OF TWO IDENTITIES 


Donald A. Steinberg 


ABSTRACT. This paper shows, by means of two examples, the possi- 
bilities of deriving expressions by the use of the notions of permutation 
and combination. One of the examples is the binomial theorem. 


By the expansion of various polynomial expressions it is possible to 
derive many identities in integers, some having practical applications, 
and others exemplifying curious properties of numbers. This method is 
utilized in Netto’s Lehrbuch der Combinatorik to obtain a large number of 
such identities. Many of these may be derived using different techniques. 
In this paper I propose to use the notions of permutation and combination 
to derive a familiar identity, the binomial theorem, and also to derive a 
new identity which, like the binomial theorem, is an expansion of a quan- 
tity raised to a power. 

1. We first see how the notions of permutation and combination can be 
used to give a simple proof for the binomial theorem. Consider a sequence 
2,2 


-,2, of n variable terms and a domain {a,,---,a,,} of m distinct 


ot 
values such that each term of the sequence can assume any of the values 
of the domain. Clearly the total number of actual sequences (i.e., an ac- 
tual sequence being obtained by giving each z,; a value from {a,,---,a@,,}) 
ism". We can also speak of this as the total number of arrangements of 
the sequence. We now partition {a,,---,a@,,} into two subsets A and its 
compliment A’, and specify that neither subset be empty. Let the number 
of elements in A be p and the number in A’ be (m-p). Let us consider the 
number of arrangements in which values from A are assumed by any & 
terms of the sequence. The expression for this is easily seen to be: 


(,)F* 


The remaining (m-p) values are assumed by the remaining (n-&) terms in 
(m-p)"-* ways. Thus the total number of ways in which any & terms can 


()e* (m—p)* 


207 


assume p values is: 





208 MATHEMATICS MAGAZINE (March-April 


If we let k=0,1,-+--,n every possible arrangement is counted at least 
once. For, every arrangement can be characterized by the number and 
positions of terms the values from A occupy. No arrangement is counted 
more than once, for that would mean values from A occupy & and k’ terms 
in the same arrangement, where & + k’. Hence: 


S (7 )e (m-p)"-* = m”. 
k=0 


Putting m-p = g, m = p+q, we have: 


n 


a (")e" q”* = (q+p)”, 


k=0 
which is the desired result. 


2. We now employ a different combinatorial situation to derive a new 
identity. Again we consider a sequence of length n and a domain of m 
distinct values. We distinguish two cases, that for which m<n and that 
for which m>n. 

First let m<n. We require that the first kterms of the sequence assume 
distinct values and that the (k+1) term assume a value already assumed 
by one of the first & terms. The remaining (n—k-1) terms have no restric- 
tions. The number of such arrangements is seen to be: 


m! km"™—*-! 


(m—k)! 


Let k=1,2,---,m. Then every possible arrangement will be counted 
once and only once. Every arrangement will be counted at least once be- 
cause every arrangement has the property that for some &, its first & and 
not more than & terms are distinct. Conversely, no arrangement will be 
counted more than once; for, no arrangement has the property that for & 
and k’, its first k and not more than k, and its first k’ and not more than 
k’ are distinct where k+k’. The number m is the proper upper limit of the 
summation. For, since the domain consists of m values we can have ar- 
rangements in which there are at most m distinct values. Thus we have: 


m m! km™—~*- 1 


me og” 
iS (m—k)! 
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m 
¥ 3 k 1 
(m—k)! m*+ 1* m! 
k=1 


If m>n the preceding discussion remains valid except for that which is 
concerned with the upper limit. In this case n becomes the upper limit be- 
cause, since there are n terms in the sequence, there can be at most n 
distinct values. Further, we note that when k=n, the number of arrange- 
ments is not: 

nm! 


m(m=—n)!’ 
which is the result of putting k=n in the general expression, but is: 


m! 
(m—n)! 


Our identity now is: 


n—k—1 ' 


n—-1 

S km! m m! - 
+ =m, 

hawt (m-k)! (m—n)! 

k=1 





I should like to express my appreciation to Dr. D. H. Lehmer of the 
University of California for his kind encouragement and guidance in the 
preparation of this paper. The publication of this paper was sponsored by 
the U. S. Atomic Energy Commission. 
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NOTE ON THE DERIVATIVES OF THE LEGENDRE 
POLYNOMIALS 


Gideon Peyser 


As is well known, see [1], the Legendre Polynomials P,(z) in the in- 
terval -1<2<¢1 assume their maximum absolute value at z=1 i.e. {P(x)| < 
P,(1) = 1. Furthermore for n>0 and z interior to the interval, |P,(z)| <P,(1). 
We shall show that these properties also hold for all the derivatives of 


the Legendre Polynomials. 


Theorem m 

1™P (zx) 
If PL™%(g) = a 
dz 


then for O<m<n 
[P™(2), < PQ) -1 sz) 


Remark : Actually, it is known, see [2], that 
Ppo™ (1) = (nem)! 
44 (n—m)! m! 2™ 


for all m and n. 
Proof 
For n=0 the theorem is correct, as P“(x) = P®(1) =1. Suppose the theo- 


rem is correct for n—1 and all m<n-1 and we shall prove it for n, For 
m=Qit is the above statement about the Legendre Polynomials. Using the 


recursion formula, see [2], 
(m) (m— 1) (m) 
Pr (a) = (n+ m—-1)PLEy (a) + @P Oh (2) 


we have for 0<m<n 


|PL™(2)| < (n+ m= VPI Ma)| + [2PM (@)| 


< (n+ m-V)PO™ (1) + P™ (1) = PO™ (1) 


For m=n we have P<"(z) = P{™(1) which concludes the proof. 
Furthermore it follows from (2) that for m<n the inequality (1) holds as 

a strict inequality for any point z inside the interval -1 < 2 <1 i.e. 

|P™(e)| < P™(1), 
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A METHOD FOR FINDING THE REAL ROOTS OF CUBIC 
EQUATIONS BY USING THE SLIDE RULE 


Louis L. Pennisi 


Methods for solving cubic equations by means of the slide rule have 
been investigated by various authors, [1] and [2]. This note will give a 
somewhat different treatment including criteria concerning the nature of 
the roots which are obtained in a simple manner from the coefficients of 


the equation. 
Suppose that we wish to employ the slide rule to solve the cubic equa- 


tion 

(1) y> + Ay? +By+C=0. 
Substituting y = z-A/3 into equation (1), we obtain 
(2) 2° 4 9qr+p=0, 


where 


(3) q=(B-1/3A") and p= (2/27A*-1/3AB+C). 


If we denote by &,, R,, R, and r,, r,, r, the roots of equations (1) and 
(2) respectively, we then have 


(4) R,=7,-4/3, R,=7,-A/3 and R,=7,-A4/3. 


Thus, our task is to first find the real roots of equation (2) by using the 
slide rule, and then use equations (4) to find the roots of equation (1). We 
shall solve equation (2) by considering the following cases. 

CASE 1. When qg is a positive number, c, and p is a negative number, 
-d, in the special cubic equation (2). 

For this case, we may write equation (2) as 


(5) a(r*+c)=d, (c>0, d>0). 


Cubic equations of the form (5) have only one positive real root, say, r,. 
The following theorem gives us useful information concerning the nature 
of this root. 
THEOREM 1. 
(i) Ifd¢l, thenr,<1. Also r.<¥a ‘ 
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(ii) If es 1, and d>2, thenr,>1. Alsor,<¥d. 
(iii) If d>c+1, thenr,>1. Also, r,<¥d. 
(iv) /f c>1 and d<c+1, thenr,<1. Alsor,< Ya. 
Proof. (i) Since r, >0 is a root of equation (5), we have r,(r? +ce)=d< I 
Since c>0, then r,<1, for, if r,>1, then r (rn? 4 c)>1, a contradiction. 
Since r,c>0, then r (rts c)=d implies that r?<d. Hence r.<¥d. 


(ii) Since d>2, then r (r5+c)>2. Since 0<c<1 and r,>0, then r,>1, 


2 
1 
< 


for, if r,<1, then r (ris ce) <2, a contradiction. The proof that r,<¥d in 


(ii), (111), (iv) is the same as given in (i). 
(iil) Since d>c+1, then rts c)>c+1. Hence, (r,-1)(r*+r,4 
Since r,>0 and c>0, then r,>1, for, if r,<1, (r,-D(ri +r, +04 


contradiction. 


2 


(iv) Since d<c+1, then r (rt+c): c+1, Hence, (r,-1)(r5+7r,+0e+1)<0. 


Since r,>0 and c>0,thenr,<1, for, if r,>1, then (r -UW(rt+r, 4 e+1)>0, 


a contradiction. 

For this case, we use the following procedure to find the only positive 
real root on the slide rule. 

PROCEDURE 1. 

(1) Move the hairline to the number d on the D scale. 

(2) Move the slide until the value of z on the C scale is under the hair- 
line and the left (or right) index of the C scale is on the value of (x?+ c) 


on the D scale. 
(3) The real positive root, x=r,, is read under the hairline on the C 


scale, 

REMARK 1. When g is a positive number c, and p is a positive number 
d,in equation (2), we may also use CASE 1 by simply finding the real root 
of f(-x) =0, where f(z) =2°+cx-d. 

CASE 2. When gq is a negative number, -c, and pis a negative number, 
-d, in the special cubic equation (2). 

For this case, we may write equation (2) as 
(6) a(zr*-c) =d, (e>v, d>0). 

Cubic equations of the form (6) have always one real positive root, 
say, r,. The following two essential theorems can be easily proven. 

THEOREM 2. /fr, is the real positive root of equation (6), then r, and 
r, are real negative roots of the same equation if and only if r?< 4/3. 

THEOREM 3. (i) /f 7, is the real positive root of equation (6), then 


riz7vec. 
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(ii) If r, is the real positive root of equation (6) and r*<4/3c, then 
Ira|<Ve and |r,|<yc. 

REMARK 2. For this case, we may check our roots by utilizing the 
following facts: 

(i) The sum of the roots is equal to 0, that is, r,+r,+7r,=0. 

(ii) The product of the roots is equal to d, that is, r,r.r,=d. 

We shall use the following procedure to find the real positive root of 
equation (6). 

PROCEDURE 2. 

(1) Move the hairline to the number d on the C scaie. 

(2) Move the slide until the value of z on the C scale is under the hair- 
line and the left (or right) index of the C scale is on the value of (x?-c) 
on the D scale. 

(3) The real positive root, x=r,, is read under the hairline on the C 
scale. 

Since the expression (x?-c) is not altered by putting x equal to -z, we 
use Procedure 2 also to find the real negative roots of equation (6). 

REMARK 3. When qis a negative number, -c, and p‘is a positive number, 
d, in equation (2), we may also use CASE 2 by simply finding the roots of 
f(-x) = 0, where f(z) = 2°-cxr-d. 

EXAMPLE. Solve y*+3y?-4.75y-10.5 = 0. 

Using (3), we find g=-7.75 and p=-3.75. Hence, (2) becomes x°~7.752 
—3.75=0 or a(z ?-7.75) =3.75, which belongs to CASE 2. By Theorem 3, 
r,>vyi.?5=2.78. From Procedure 2, we have: Move the hairline to 3.75 
on the D scale. Move the slide until 3 on the C scale is under the hair- 
line and the left index of the C scale is on 1.25 on the D scale. Hence 
r,=3, since when z=3, 2?—7.75=1.25 which is the same value found at 
the left index of the C scale on the D scale. From Theorem 2, since 9< 


4/3 (7.75) =31/3, r, and r, are real negative roots. From Theorem 3(ii), 
|r,|<2.78 and |r,|<2.78. Always keeping the hairline on 3.75 on the D 


scale, we use Procedure 2 twice to find r, and r,. Move the slide until 


2.5 on the C scale is under the hairline and the left index on the C scale 


is on 1.5 on the D scale. Hence r,=-2.5. Similarly we find r,=-0.5. From 
Remark 2, we see that 3-2.5-0.5=0 and 3(-2.5)(-0.5)=3.75. Finally, 
using (4) we find R ,=2, R,=-3.5 and R,=-1.5. 
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LINEAR DIOPHANTINE EQUATIONS 


Arthur B. Brown 


We shall describe elementary methods for solving linear Diophantine 
equations, either single equations or systems of equations. 

The equations with which we shall deal] have integral coefficients and 
require integral solutions. Such equations are named after Diophantus of 
Alexandria, who lived probably between 100 and 300 A. D. 

The methods are explained with the help of specific examples. We 
begin with the case of a single equation. 

Method I If an equation has an unknown with coefficient 1 or —1, then 
all the other unknowns can be taken as arbitrary integers, and the parti- 
cular unknown is then determined from the given equation. 

Example 1. 4e+7y-2+3w = 18. The solution is: Take z, y, w as ar- 
bitrary integers, and determine z from the given equation. It is clear that 
all solutions in integers are thus obtained. 

Method II If the numerically smallest coefficient, say a, is notlor 
+1, but is a factor of each of the remaining coefficients, and is not a fac- 
tor of the right member of the equation (the “constant”term), the equation 
is incompatible, that is, it has no solution in integers. If a divides all the 
coefficients and also the right member, we divide through by a, thus ob- 
taining an equivalent equation to which Method | can be applied. 

Example 2.8x+14y—-22+6w = 36. If we divide by 2 we get the equation 
of Example 1 and can apply Method I. If the right hand side were 37 in- 
stead of 36, the equation would have no solution in integers, as 37 is not 
divisible by 2. 

Method II can also be applied to any single equation of a system of 
simultaneous equations. If one equation is thereby shown to be incom- 
patible in integers, then the system is incompatible in integers. 

We now continue with the case of a single equation. 

Method Ill. If the numerically smallest coefficient, say a, is not a 
factor of all of the remaining coefficients, the equation can be replaced by 
a new equation with a non-zero coefficient numerically smaller than a, 
by the method used in the following example. 

Example 3. 164+ 42y-21z2 = 100. 

We first set 


r= 2°-3y+2. 
215 
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The coefficients on the right in (2) are determined as follows. The coef- 
ficient of z’ is 1, which is the rule for the case that z has the numerically 
smallest coefficient in (1). The coefficients -3 and 1 of y and 2 respec- 
tively in (2) were chosen so as to minimize the numerical values of the 
coefficients of y and z in the equation obtained by substituting from (2) 
into (1). That equation is 


l6z’-6y-52 = 100. (3) 


After solving (3) in integers, we shall use (2) to determine z. 
Applying Method III to (3), we set 


B= 2°-y+32’, 


which converts (3) to 


, 


2’-y-5z’= 100 
which has solution 
x’, 2” arbitrary integers, 
y = z’-52°-100. (7) 
Substituting into (4) we get 
2 = 22°+62°+100. (8) 


Then, from (2), (7) and (8) we have 
z = 212+ 400. (9) 


All solutions of (1) in integers are given by (6), (7), (8) and (9). 

Application of Method III, repeated as often as necessary, must even- 
tually lead to an equation all of whose coefficients are integral multiples 
of a particular coefficient a. We then apply Method I or Method II, depend- 
ing on whether or not a = 1. 

For an equation in two unknowns, Method III, combined with I and II, 
is roughly equivalent to the standard treatment of a Diophantine equation 
in two unknowns.* When the number of unknowns is greater than 2, substi- 
tutions such as (2) above bring the problem to solution more quickly than 
the standard methods. 

We turn now to systems of linear Diophantine equations, giving pro- 
cedures simpler and more efficient than familiar ones. 

Method IV .If a coefficient 1 or—1 appears in an equation of the system, 
the system can be replaced by a system in one less equation and in one 
less unknown. (The number of unknowns that actually appear with non-zero 


*See Uspensky and Heaslitt, Flementary Number Theory, New York 1939, pp. 52-64. 
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coefficients may, however, decrease by more than unity.) We illustrate 


with the system 


32+ 2y-52+ 4w = * (10) 
4e- y+T72 (11) 
5a + 6y -2w = 9. (12) 


Noting coefficient -1 for y in (11), we add twice the left member of (11) 
to that of (10), and similarly on the right, obtaining 


llz+92+4w = 24. (13) 
Similarly, from (11) and (12) we get 
292 + 422-2w = 69, (14) 


The given system of 3 equations in 4 unknowns is now seen to be equiv- 
alent to (11), (13), (14). But as (11) can be used to determine y when z, z 
and w have integral values satisfying (13) and (14), we thus see that the 
problem is essentially reduced to that of solving the system (13), (14), a 
system of 2 equations in 3 unknowns. 

Method V. If no coefficient is 1 or ~1, but the coefficients of some un- 
known are relatively prime, the system can be replaced by one having a 
coefficient 1 or -1. We illustrate with the system 


62+ ++. = (15) 
152 + ++. (16) 
202 + +++ ° (17) 


By adding an integral multiple of each side of an equation having nu- 
merically smallest coefficient, here 6, of the given unknown, here z, to the 
corresponding side of one of the other equations, we get an equation with 
coefficient of z numerically less than 6. Thus, multiplying (15) through 
by -2 and adding to (16) we get 


324+--- = b-2a. (18) 


We could similarly work with (15) and (17), but we note that (15), (16), 
(17) are equivalent to (15), (17), (18), and by multiplying (18) through by 
-7 and adding to (17) we get an equation with coefficient -1, namely 


—%+--. = c-75+4 14a. (19) 


The given system (15), (16), (17) is equivalent to (15), (18), (19). 
If a system has no coefficient 1 or -1, and the coefficients of each un- 
known have some common factor, while the coefficients in each equation 


have no common factor, then none of Methods II, IV, V can be applied. In 
such a case we use 
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Method VI. Under the circumstances described in the preceding para- 
graph, we apply Method III to an equation of the system having a non-zero 
coefficient numerically as small as possible, but we make the substitution, 
such as (2), in all the equations of the system. For example, given the 
system 


4r+ 9y+10z2- 2w (20) 
20r+ 6y- 52+14w=7 (21) 
6r+18y+ 52- 4w = 10 (22 


we note that -2is the numerically smallest coefficient and accordingly 
make the following substitution suggested by examination of equation (20)* 


w= w’+4y. (23) 
This results in the equations 


Avr+ y+10z2- 2w’=9 (24) 
202+62y- 52+14w’=7 (25) 
6r+ 2y+ 52- 4w’= 10, (26) 


a system having a coefficient unity. The problem of solving (20), (21), 
(22), is now reduced to the problem of solving (24), (25), (26), since if z, 
y, 2, w’ satisfy the latter system, w can be obtained from (23), and (20), 
(21), (22) will be satisfied. All solutions in integers are obtained, since 
if 2, y, 2, w satisfy (20), (21), (22) then 2, y, z, w’ satisfy (24), (25), (26). 

In general, after getting a new system such as (24), (25), (26), we apply 
Method II, IV or V if possible (and here of course Method IV would be 
used next); otherwise continue with Method VI. 

It is clear that use of Methods I to VI will result in obtaining all solu- 
tions in integers of the given system of equations, and reveal incompati- 
bility, if there is such, when an equation is reached having a common fac- 
tor of the coefficients on the left which is not a factor of the number on 
the right side of the equation. 

Examination of Methods I through VI shows that arbitrary integers are 
introduced only when the work is at a stage where a single equation is to 
be solved, and when Method I is applied. The number of arbitrary integers 
is one less than the total number of unknowns that must satisfy the single 
equation if not all the coefficients are zero, and is the same as the num- 
ber of unknowns if all the coefficients are zero (that is, the equation in 
question is “0=0”). In either case, the number of arbitrary integers equals 
the number of unknowns minus the rank of the matrix of coefficients, for 
the single equation. Now, examination of Methods II through VI shows 


*We might also use w=w’+2r+4y+5z, thus replacing (20) by a much simpler 
equation, but (21) and (22) would not be simplified. 
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that in each case an equation, or a system of equations, is replaced by an 
equation or system for which the quantity “number of unknowns minus rank 
of the matrix of coefficients” is unchanged. Hence the following theorem, 
except possibly the final sentence, is valid. 

THEOREM. /f a linear equation or a system of linear equations with 
integral coefficients is compatible in integers, all solutions in integers are 
given by equating the unknowns to linear functions, with integral coef- 
ficients, of n-r arbitrary integers, where nis the number of unknowns and 
ris the rank of the matrix of coefficients. No fewer than n—r arbitrary in- 
tegers will suffice. 

We now prove the final statement of the theorem. Let (S) denote the 
given system of equations and ¢,, ¢,, +++, ¢, the set of arbitrary integers 
on which the solutions in integers of (S) depend, so that & = n-r. Suppose 
that all solutions in integers of (S) could be given by 


Bp = Cy My t+ Cg Mater + Cm Uy td; (j= 1,2,+++,n) (27) 


with u,, U,, *++, U,», integers, and m<k. Now the formulas for @ i, es, #2, in 
terms of ¢,, +++, ¢, are the result of a sequence of formulas similar in na- 
ture to (2), that is, the new unknown is always introduced on the right 
with coefficient unity. We consider (2) to give formulas for z, y, 2 in terms 
of x2’, y, 2, so that the equations y=¥, z2=z2 may be introduced to complete 
the set of formulas for the change. Thus (2) with the two additional equa- 
tions can be solved, giving 2’=2+3y-z2, y=y, z=2, that is, the new un- 
knowns in terms of the old unknowns. A similar statement can be made in 
the general case. Since a similar situation holds at each introduction of a 
new set of unknowns (where of course n-1 of the equations are simply of 
the form “unknown equals itself”), we infer that the final set of unknowns 
introduced are linear functions of z,, +++, z,. In particular, ¢,, ¢,, «+, ¢, 
are linear functions of z,, «++, zp, 


tea fy@ytortly ey (j= 1,2,-++,%) (28) 


a fact which is not altered by the relationships which the final set of un- 
knowns (namely ¢,, «++, ¢, and n-k other unknowns) must satisfy in order 
to make z,, -+-, x, satisfy (S). From (27) and (28) we get 


bp = Jyh yte t+ Im Mth; (j= 1,-++,k) (29) 


where the g’s and h’s are constant integers, with (29) satisfied by (¢) and 
(u) whenever both give the same solution (z) of (S) in integers. 

If the u’s are allowed to take on all real values, then (29) considered 
as parametric equations of a locus in ¢,, ---,¢, euclidean space is a hyper- 
plane of dimension at most m<k. Then all k-tuples of values of ¢ ,, +», ¢, 
corresponding to solutions of (S) in integers would be confined to this 
hyperplane, which is impossible since ¢,, ---, ¢, can be arbitrary integers 
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and hence determine points not confined to any plane of dimension less 
than & in k-space. We infer that the theorem is true. 

Conditions for compatibility of (S) can be found in a paper by Oswald 
Veblen and Philip Franklin. * 


*“On Matrices whose Elements are Integers”, Annals of Math. (2), Vol. 23 (1921), 
pp. 1-15; also as Appendix II of Veblen’s “Analysis Situs”, second edition. 
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ARITHMET ICAL CONGRUENCES WITH PRACTICAL 
APPLICATIONS 
(Popular Article ) 


Carol Law 


When we say that something is congruent, we mean that it is in agree- 
ment or in harmony with something else. In geometry we speak of con- 
gruent triangles and designate their harmonious qualities by the equal and 
similar signs =. In theory of numbers the term congruence still has the 
same connotation in that groups of numbers have some relation to each 
other. We are often concerned here with properties of numbers that are 
true for a whole group of integers. For example, we find a certain property 
holding for squares of odd integers when divided by 8 leave 1 as a re- 
mainder. This properiy holds for all odd numbers, or a class of numbers 
differing from each other by multiples of 2. Another example is the fact 
that when we square any number having 6 as the last digit, the product 
will also end in 6. Thus, we see another property holding for a whole 
group of numbers. The consideration of these properties holding for all 
integers differing from each other by a multiple of a certain integer leads 
to a notion of congruences. 

Carl Friedrich Gauss, who lived in the 18th century, is attributed with 
the foundation of congruences. He was a mathematical genius, especially 
with theory of numbers. It is this topic, arithmetical congruences, which 
he treated so well, that will concern us now. Just what do we mean by 
arithmetical congruences? 

In my paper I will attempt to do three things: (1) show what is meant 
by two numbers being congruent, (2) show that it is possible to solve for 
unknowns in linear, quadratic, and cubic congruences, and (3) show how 
we might apply congruences practically. 

What do we mean when we say that two numbers are congruent? Two 
integers a and 6 whose difference a-b is divisible by a given number m 
(not 0) are said to be congruent for the modulus m or simply modulo m. 
This means simply that two numbers are congruent if they have the same 
remainder when divided by the same number. This is expressed a= 6 (mod 
m) where the congruent sign is written =. For example, 17=5 (mod 12) or 
9?=1 (mod 5), Note that the congruence is true either if you substract the 
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integers first or reduce them first. For example, 
Substracting first: Reducing first: 


9?_1=0 (mod 5) 9?=1 (mod 5) 
81/5 = 1/5 (mod 5) 
81-1=0 (mod 5) 1=-4 (mod 5) 
0=0 (mod 5) 1=1 (mod 5) 


Notice also that -4 and +1 are the same remainder modulo 5 as can be 
shown by dividing 1 by 5 as follows: 


+ 


A geometric interpretation of congruences is helpful. The first figure 
is a representation of integers, and the second figure represents these 
integers modulo 6. 





alk -3 -2 -] 0 


Figure I 


Figure II 


Notice that the numbers at any one point in Figure II are congruent with 
each other modulo 6. For example, 


-5 = 7 (mod 6) 
-5-7 = 0 (mod 6) 
-12 = 0 (mod 6) 
0 = 0 (mod 6) 


An application in daily life of this geometric interpretation is the 
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clock. The hands indicate the hour modulo 12. Below is a representation 


of army time. Whereas we start over with 1:00, 2:00, etc., after reaching 
12:00, soldiers begin with 0100 and continue up to 2400 for an entire day. 


1500 10300 0900 } 2100 


Notice at each position 12 is added so as to make congruent numbers. 


0100 = 1300 (mod 12) 
0100- 1300 = 0 (mod 12) 
- 1200 = 0 (mod 12) 

0= 0 (mod 12) 


Many properties of ordinary equations hold for congruences. Some of 
these are as follows: 


. For any modulus m, a = a (mod m). 

. If a = b (mod m), then 5 = a (mod m), 

. If a = b (mod m) and } = c (mod m), then a = c (mod m). 

. If a = b (mod m) and c = d:(mod m), then a = c = b = d (mod m). 

. If a = 6 (mod m) and c = d (mod m), then ac = bd (mod m). 

. If a = b (mod m), then a” = 6” (mod m) where n is a positive integer. 

. A common factor n cannot be removed from na = nb (mod m) unless it is 
a prime of m or unless m is also divided by n. 

. If a = 6 (mod m) and d is a divisor of m, then a = 5 (mod @). 

. If a = b with respect to moduli m,, m,, --- m, then a = 6 (mod /)where/ 


is the least common factor of m,, m,, -+- m,. 
Let us see how we can apply these principles. Suppose we wanted to 
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find if 47 is a factor of 27*-1. To obtain 27° would constitute a good bit 
of work, but with congruences the solution is simple. We simply break 
down 27° into powers of 2 that we can readily compute, divide each term 
by 47, multiply the remainders together, again divide by the modulus, and 
end up with the final remainder. Since the remainder equals 1 we can say 
that 47 divides 27-1 evenly. This problem is shown as follows: 


= 1 (mod 47) 

= 64-647.32 (mod 47) 
17- 7 +32 (mod 47) 
3808 (mod 47) 
1 (mod 47) 
0 (mod 47) 


Arithmetical congruences become more valuable with the solution of 
linear, quadratic, and cubic congruences with one unknown. Actually, the 
solutions to each are similar. For example, the linear congruence is 
solved in the following manner: 22+1 = 0 (mod 5). The method is simply 
by trial and error—substituting numbers for z and seeing if they satisfy 
the congruence. However, there is a limit to the amount of numbers that 
need to be substituted. Mathematicians have found that only % of the 
positive and negative values of the modulus m need to be tried. Higher 
numbers, as I have already indicated, can be shown to be the same re- 
mainder as one of the lower values, In this case all values up to and in- 
cluding + 2 need to be tried. We find + 2 is a solution as well as - 3; but 
+2 and -3 is the same remainder modulo 5. 

This method holds for the solution of a quadratic also; however, if the 
congruence can be factored, the work involved is reduced. For example: 


2*+52+6 = 0 (mod 7) 
(x+3)(x+2) = 0 (mod 7) 
(2 +3) = 0 (mod 7) 

#« = 4 or -3 (mod 7) 
(2 +2) = 0 (mod 7) 


xz = 5 or -2 (mod 7) 


For a cubic congruence such as 2*+1927+2+23 = 0 (mod 42), one can 
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readily see that to substitute in 44 the positive and negative values of the 
modulus 42 would be quite monotonous, so mathematicians generally em- 
ploy the Chinese Remainder Theorem to solve the cubic. 

Usually the human mind attempts to find some practical application of 
this knowledge. I have developed an original idea which | think is inter- 
esting and perhaps useful to some, in my case to the marching-band di- 
rector. Let us take a very simple quadratic congruence and illustrate how 
it might be employed in figuring marching-band formations. 

Let us suppose that the college band of St. Benedict’s College were 
marching down Pittsburg’s football field. Mr. Roark, St. Benedict’s band 
director, decided to have his band march on the field in a perfect square 
formation, that is as many rows as columns. Then at the middle of the 
field he wanted them to go into a P formation for Pittsburg. To make a P 
he wanted 9 men for the loop and 15 for the stem. He wanted to enclose 
the P with rows of men, 15 per row. How many men would he need? We 
could set up the following relation: 


2?~9 = 0 (mod 15) 
(2-3)(2+3) = 0 (mod 15) 
(2-3) = 0 (mod 15) 

+3 or -12 = 0 (mod 15) 
(x+3) = 0 (mod 15) 

-3 or +12 = 0 (mod 15) 


Using 12 as a solution, we have, 


144-9 = 0 (mod 15) 
135 = 0 (mod 15) 


Now we have 12 rows and 12 columns of men marching down the field. 
Nine of these will be used for the loop of the P which leaves 135 for the 
stem and the rows. By dividing 135 by 15 we see that there are 9 rows of 
15 men in each. We can therefore set up the following formation: 


Perhaps one of the most practical and fascinating applications of 
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congruences is in working calendar problems. 

You have heard about people being “walking dictionaries.” You can 
become a “walking calendar.” We can figure out the days of the week on 
which certain dates fall; certain dates of the month on which certain days 
fall; the day on which you were born; the day of your next birthday; or 
even dates on which Easter falls. One day in history class when we were 
discussing the second World War, my professor asked me, “When did the 
United States declare war on Japan?” I hesitated momentarily, then an- 
swered emphatically, “On Monday, December 8, 1941.” He looked at me a 
little surprised to think that I remembered even the day of the week; how- 
ever, he really would have been surprised if he had known that I didn’t 
remember but that I quickly applied a mathematical formula. This is what 
I did: 


f = ?, day Pee ee pe ie MH, (mod 7) 
k = 8, date 44 
m= 10, month f = 8+(2.6)(10)-0.2+41 ne (ae (mod 7) 
D = 41, years in century its 
C = 19, century 
f = 8+25+41+10+4-38 (mod 7) 


f = 88-38 (mod 7) 
f= 50 (mod 7) 
f= 1 (mod 7) 


Since Sunday is 0, Monday is 1,so Monday is the required day of the week. 

We could apply congruences in many ways in the home but let’s see if 
we can’t apply them in the social world. Let’s say you’re coming home 
from the midnight ball with your fiancée. The moon is shining like a silver 
dollar—and the subject, guess what! Finally you agree—yes, June is the 
month and Saturdays are the only days you'll consider. “But what date 
shall we set? Of course you have no calendar with you and at such a vital 
moment! But, yes, she had studied congruences. She’d figure out on what 


dates Saturdays fall! 


C 
f =-1, day <k = (2.6m-0.2)-2a-46+ | -C—f (mod 7) 
m= 4, month 19 
a=1 ~k = (2.6)(4) -0.2-2(1)-4(4)+ ~ ~19-(-1) mod 7) 
4 
b=4 
C = 19, century -k = 10-2-16+4-19+1 (mod 7) 
k = ?, date -k = ~22 (mod 7) 
kel (mod 7) 
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Note: a and 6 are obtained by dividing the year by 4 and 7 respectively 
and using the remainders. (When using these formulas always use the 
least whole integer.) She found that in June, 1957, Saturdays fall on the 
first, the eighth, the fifteenth, the twenty-second, and the twenty-ninth. 
The remainder 1 equals the first date of the month. Simply add 7 dayson 
to this date each time to find the remaining Saturdays. She applied Zeller’s 
Rule and the wedding date was set! 


REFERENCES: 


Elementary Number Theory by Uspensky and Heaslet 
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PROBLEMS AND QUESTIONS 
Edited by 
Robert E. Horton, Los Angeles City College 


Readers of this department are invited to submit for solution problems believed to be 
new and subject matter questions that may arise in study, in research, or in extra-academic 
situations. Proposals should be accompanied by solutions, when available, and by any 
information that will assist the editor. Ordinarily, problems in well-known textbooks 


should not be submitted. 
Solutions should be submitted on separate, signed sheets. Figures should be drawn 


in India ink and twice the size desired for reproduction. 
Send all communications for this department to Robert E. Horton, Los Angeles City 


College, 855 North Vermont Ave., Los Angeles 29, California. 


PROPOSALS 


334. Proposed by Huseyin Demir, Kandilli, Eregli, Kdz, Turkey. 
Find the simplest expression for the area S enclosed by the arc AM of 
a cycloid, the arc TM of the rolling circle 2 (a) and the base line segment 


AT. 


335. Proposed by Robert E. Shafer, University of California Radiation 
Laboratory. 
Prove 


n 


5 sat) > 


t=] t=] 
for all n > 1. 


336. Proposed by C.W. Trigg, Los Angeles City College. 

A uniform bar with rounded ends and length z is hung from one end by 
a string of length z and negligible mass. The other end of the string is 
attached to a vertical wall. When the free end of the bar is placed against 
the wall, it is found that @ is the smallest angle that can be made between 
the bar and the wall so that the bar will not slipand fall down. (The plane 
of the bar and string is perpendicular to the wall.) What is the coefficient 
of friction between the bar and the wall? 


337. Proposed by Victor Thebault, Tennie, Sarthe, France. 
Determine the right triangles whose sides are integers and whose hy- 


potenuse is twice a square. 
228 
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338. Proposed by M. S. Klamkin, AVCO, Lawrence, Massachusetts. 


R R R 
The vector field > satisfies the equations V *-> = 0 and V -~> =0. 
r r 


. r 
Consequently, this field has a scalar potential and a vector potential. The 


1 
scalar potential is well known to be —.Determine the vector potential. 
r 


339. Proposed by D.F.. Huntington and D.A. Kearns, University of Maine. 
Criticize the following “proof” that the Cantor Ternary Set (Cantor Dis- 
continuum) is non-denumerable. 
At the nth stage of the construction of the set there are 2” closed in- 
tervals, the 2"*' end points of which belong to the set. Since there are a 
denumerable number of stages in the construction, the total number of end 


‘ sett. : 
points is 2 = 2 “ =5N,. Therefore these end points are non-denu- 
merable, and since they constitute a subset of the Cantor set, the entire 
Cantor set is non-denumerable. 


340. Proposed by R. G. Buschman, University of Wichita. 
Prove that 


¢(k) = af ({e’] —e° + 1/2) dt 
0 


— (-1)*(k-1) 
e k(k+1) 


where [z] represents the greatest integer not exceeding z and ¢(k) is the 
Reimann Zeta-function. 


SOLUTIONS 


Christmas Sunday 


313 [September 1957] Proposed by Sidney Kravitz, East Paterson, New 
Jersey. 

Show that Christmas falls on a Sunday more often than once every 
seven years. 

Solution by C.W. Trigg, Los Angeles City College. A usual year con- 
tains 365 or 7(52) + 1 days, and every year divisible by 4 contains an 
extra day unless the year is divisible by 100 and not by 400. It follows 
that 400 years will include 97 leap years and exactly 146097 days or 
20871 weeks. Hence, the calendar repeats every 400 years. 

Christmas will fall on Sunday in 1960 and will progress through the 
seven-day week one day at a time for three days and then will skip a day 
in leap year. The pattern will repeat every 28 years until broken up by a 
multiple of 100. For example: 
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Mon. Tues. od. Thurs. Fri. 


1961 1962 

1967 1969 

1972 1973 1975 
1879 1980 


2302 2303 


Now Christmas falls on Sunday only once in the span 1961 to 1967, 
inclusive. Hence, it does not fall on Sunday more often than once “every 
seven years”. 

If, ‘however, the statement is interpreted to mean that “the probability 
of Christmas falling on Sunday is greater than 1/7”, then we merely need 
to consider the days upon which Christmas falls in the 1960-2359 span. 


Sun. Mon. Tues. Wed. Thurs. Fri. Sat. 
58 56 58 57 57 58 56 


Considered in this light the proposition holds. 

A -similar problem in AMERICAN MATHEMATICAL MONTHLY, 40, 
607, ‘December 1933 shows that the thirteenth is more likely to fall on 
Friday ‘than on any other day of the week. 


Also solved by George Bergman, Stuyvesant High School, New York 
City; ‘Sam Kravitz, East Cleveland, Ohio; and the proposer, all using the 
probability approach. D. A. Breault, Sylvania Electric Products, Inc. 
Waltham, Massachusetts, showed the proposition false between 1967 and 
1976. 


Powers of Two 


314. [September 1957] Proposed by J.M. Gandhi, Lingraj College, Belgaum, 
India. 

Show that 2"-' = 2(,C,-,C,) + 4,C,-,C,) + 6(,C,-,C,) + +++; the 
last termbeing n(,,C,,) when n is even and -(n-1),,C,, when n is odd. 


1. Solution by W.B. Carver, Cornell University. Consider the case when 
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n is odd. We write C, for ,C, and /, for ,_,C,. It is well known that 
27-1! = 2/7 + 2/, + 2/, + 2/7 «+ + 2/7_,. We may write this 


2 1-2/5 +(4-2)/5+(6-4)/; «++ + {(n-1)-(n-3)}/)_ 
=2(/5-15)44(/q-1¢)46(/ 5-17) -+4(0-3) Sp g/g) +(0-D)_n 
=27 +1 9-1-1 gt 4/411 15) 00 10-3) gt g-gn g-/ 2) + 
(n—1)/y 9» 
=2(C ,-C ,)+4(C ,-C ,)+6(C .-C ,) +++ + (n-3)(C,,_-Cy_,) + (n-D)/,_,- 
Since ,_, = C,_,-C,, this proves the theorem for n odd. A slight modifi- 
cation of the above takes care of the case n even. 


Il. Solution by Michael J. Pascual, Burbank, California. Consider the 


function 
(1-2?)” 


(1) flw)=-1/r+ 
fle) =-1/2+1/2-,C + ,C 2° -,Cyr>+---+(-1)",C 27! 


f(a) =-,C ,+3,C,@7-5,C 24 +--+ (-1)"(2n-1), 0.27"? 


f(a)=- ,C ,+3,C,-5,C, +++ + (-1)"(2n-1),C,, 


and by differentiating (1) we easily have that /’(1)=1=,,C, hence 


Co =~ nl 1 +3,0,-5,C +e + (-1)"(2n-1),,C,, 


n 


C+ stn’ ott nO n=4yF 9-4 € 9+8,,C 8,6 gt-+-+L(-1)"(2n-1)+1] me 


n-on 


but the left member is (1+1)” = 2”. Hence 


-1)"(2 
BP 4 oe HO eC HMC LC ecclas ee nc 


n-2n~ 8 n-4n-s 2 





(-1)"(2n-1)+1 
> =n 





and if n is even 


oo (-1)™(2n-1)+1 ; 
whereas if n is odd > = —(n—1) as required. 





Also solved by Harry D. Ruderman, New York, New York; Huseyin 
Demir, Kandilli, Eregli, Kdz., Turkey; E. P. Starke, Rutgers University; 
C.W. Trigg, Los Angeles City College; Chih-yi Wang, University of Minn- 
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esota and the proposer. 
Equivalent Integrals 


315. [September 1957] Proposed by P. D. Thomas, Eglin Air Force Base, 
Florida. i. 
Under the restriction f(c-—z) = f(x) show that Kc bx) f(x) dz may be 


written (a+ (80)/2)( fa) dz. 


I. Solution by Ben K. Gold, Los Angeles City College. Since fic-2)= 
f(z) we have 


c c 
[ilie-a) + cl f(c-2)[-d(c-zx)] = -[2fedde + of fladde 


[pare = sho f(x) dz 


fo lovte) fladde = (a+be/2){ fla)de 


This leads to 


The result 


follows. 


Il. Solution by Waleed A. Al-Salam, Duke University. Since f(c-z) = f(z), 


2 
we have [) flalde = afr" fla) dz 


Now 


{, (a+bz) flx)dz -(" (a+bz) flx)dz a (a+be ) f(x) d. 
0 “Jo c/2 “ 


In the second integral of the right hand side, put 2 = c-y. We get 
c c/2 e/2 
Lo (a+ba) f(x) dx {9 (a+by) fly) dy ah (a+b(c-y) fly) dy 
c/2 
- 2arb0){ f(x) dz 


be 
= (a+ =f fix) dx 


Also solved by W. B. Carver, Cornell University; Huseyin Demir, 
Kandilli, Eregli, Kdz., Turkey; J. M. Gandhi, Jain Engineering College, 
Panchkoola, India; Donald A. Kearns, University of Maine; Joseph D. E. 
Konhauser, Haller, Raymond and Brown, Inc., State College Pennsylvania; 
M. Morduchow, Polytechnic Institute of Brooklyn; F.D. Parker, University 
of Alaska; Michael J. Pascual, Burbank, California; Arne Pleyet, Trol- 
hattan, Sweden; A. K. Rajagopal, Indian Institute of Science, Bangalore, 
India; Lawrence A. Ringenberg, Eastern Illinois State College; Harry D. 
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Ruderman, New York, New York; William F. Russell, St. Joseph’s College, 
Pennsylvania; Chih-yi Wang, University of Minnesota and the proposer. 


Factorization of 1 + cosné 


316. [September 1957] Proposed by A. K. Rajagopal, Lingraj College, 
Belgaum, India. 

Prove that (1 + cosn@), n an integer, has a factor (1 + cos @) if and 
only if n is odd. 

I. Solution by Joseph D E. Konhauser, State College, Pennsylvania. 
It is well-known that 1 + cos 7n@ can be written as a polynomial of degree 
n in cos 6. Let z = cos 6, and denote the polynomial by p(x). The problem 
is to show that 1 + 2 is a factor of p(x). This is the case if and only if -1 
is a zero of p(x), which is equivalent to 0=7 being a zero of 1 + cosné. 
But 0=7 is a zero of 1 + cosné@ if and only if n is odd, proving the de- 
sired result. 


Il. Solution by George Bergman, student at Stuyvesant High School, 
New York, New York. Let z = e9. Then (1 + cos 9) and (1 + cosn@) are 
equal respectively to 4(2+2+2~') and H2"+2+27"). Now an algebraic 
factor of such an expression is equivalent to a trigonometric factor if and 
only if it takes on real values for any z on the unit circle. Since this is 
the case, we need only test as to whether the first of our algebraic ex- 
pressions is a factor of the second. The factorization of the first is Ha+ 
1)(z~'+1) which the Factor Theorem shows us divides the second only 
for odd n. 


Also solved by Waleed A. Al-Salam, Duke University. Walter B. Carver, 
Cornell University; Huseyin Demir, Kandilli, Eregli, Kdz., Turkey; M. 
Morduchow, Polytechnic Institute of Brooklyn; F. D. Parker, University of 
Alaska; Michael J. Pascual, Burbank, California; Arne Pleyet, Trolhattan, 
Sweden; Lawrence A. Ringenberg, Eastern Illinois University; Harry D. 
Ruderman, New York, New York; Dmitri Thoro, University of Florida; P. 
D. Thomas, U.S. Coast and Geodetic Survey, Washington, D. C.; Chih-yi 
Wang, University of Minnesota and the proposer. 


An Exponential Inequality 


317. [September 1957] Proposed by Ben K.Gold, Los Angeles sis Ccllege. 


Prove that (e+ 2)°°-™ > (e-2)°°*™ forO0<a2<e. 
I. Solution by E. P. Starke, Rutgers University. Let 


f(x) = (e-2) In(e+z) - (e +2) In(e-2). 


Then 
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f(z) -().(—)- In(e2—2?) 

e-z e+2 
Now, remembering that .n+1/n>2 for all positive n41, it is easily seen 
that f(z)>0 for 0<|z|<e. Hence f(z) is an increasing function over the 
same range. Thus, since f(0) = 0, (e-2) Inle+z)>(e+2) Infe-2), 0<a2<e, 
and this is equivalent to the required result. 


Il. Solution by the proposer, To prove (e+zx)‘°~*>(e-a)°+™ we have 
the equivalent statement 





(e+2)° | (e-2)* a (Sct) »(222} 


e-2z e-2 


(e+ax)” (e-z)” 


which is of the form N® >N* or e>z. But our hypothesis is that 0<z2<e so 
the original statement holds. 


Also solved by J. M. Gandhi, Jain Engineering College, Panchkoola, 
India; Joseph D. E. Konhauser, State College, Pennsylvania; Michael J. 
Pascual, Burbank, California; A. kK. Rajagopal, Indian Institute of Science, 
Bangalore, India; Harry D. Ruderman, New York, New York; Chih-yi Wang, 
University of Minnesota. 


Limit of a Quotient 


318. [September 1957] Proposed by Chih-yi Wang, University of Minnesota. 
gics x 
Evaluate lim 


@0o0 (log x)* 
|. Solution by D.A. Breault, Sylvania Electric Products, Inc., Waltham, 


Massachusetts. Let 
(z)!oe% 
F(z) = l 
. (log x)” 


and we wish to find lim F(z). Set 2 = e“, and (1) becomes, 


Loco 


ew 
lim G(u) = lim =e 
Bear -nae” 
Taking the natural logarithm of both sides of (2), we see that 
lim log G(u) = lim {u? - e“log u} (3) 
U-co U0 


which clearly is -«. From this fact, we deduce that lim G(u)=0, and hence 


lim F(z) = 0, also. U-roe 
U-roceo 
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Il. Solution by M. Morduchow, Polytechnic Institute of Brooklyn. L’Hopi- 
tal’s Rule will be applied. First the derivative of the numerator and of the 


denominator must be obtained. To differentiate y=2'°*”, it is noted that 

log y=(logz)?, and y’/y=2(logz)/z; hence y’=y[2(logz)/z]. Similarly, if 

2=(logz)”, then one finds 2’= a[log (logr)+(1/logz)]. Thus Jim 2= lim -U, 
i . Hence lim ¥(1-U)=0. But U+0 as 2+ 

log (log x) + (1/log 2) £00 Z 

[as can be readily seen by noting that (logr)/z+0]. Hence y/z+0 as 2+, 





where U= 


Also solved by Walter B. Carver, Cornell University; Joseph D. E. 
Konhauser, State College, Pennsylvania; Michael J. Pascual, Burbank, 
California; Harry D. Ruderman, New York, New York;E.P. Starke, Rutgers 
University and the proposer. 


A Trigonometric Identity 


319. [September 1957] Proposed by Barney Bissinger, Lebanon Valley 


College, Pennsylvania. 
sinn6é Cos @ Cos 26 Cos (n-1)0 
a 


= cos™ ' 641 one 
*Cos@ Cos*6 * Cos 6 
values of @ for which the terms are defined. 


for those 





Show that 


Il. Solution by Lawrence A. Ringenberg, Eastern Illinois University. 
sin 20 


The equation 


cos @ 
=Ccos (1 + ) follows immediately from the double 
cos 6 


angle formula sin2@=2sin6cos 6. Suppose & is a positive integer such 


sin ké k-1 k-2 
that =Ccos 6 + cos 6cos 6 + «+. + cos 8cos (k-2)6 + cos (k-1) 6. 


sin 
Then 
sin(k+1)@ sin(k+1)@ sinké 
sind z sink@ sind 
cos k@ sin @\sin ké 
= (cos 0+ . ; 
sinké / sin 6 








=cos 6 (cos*- 19, cos*~?6 cos 6+-+-4c0s (k-1) 6) + cos ké 





cos 0 cos ké 
+o OFS 
cos 0 cos” @ 


The desired result follows by mathematical induction. 


= cos* a(1 + 


Il. Solution by Dmitri Thoro, University of Florida. Let u(x)=cosz0sec” 6, 
From the calculus of finite differences it may be readily verified (by dif- 


ferencing) that A~* u(x) = cot @sinz@ sec” 6. Hence 
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n—1 n 
S u(x) = A~* uz) = cot Osinn@sec" 6, 
2-0 z=0 


sinné | 
ca” '* @il+ 


ns Par 2 
sind [ eos cos? @ cos 


cos if cos 20 cos (n—1)0 
nl 9 


Ill. Solution by Arne Pleyet, Trolhattan, Sweden. If R(x) indicates the 
real part of z then we have 
cos@ cos(n-1)6 ‘0 et(n—1)0 | 


e 
ae n—1 aE a 
=Rcos oh feed 51g 


ooo+ 
cos. cos“ Dg cosé 


ate cost of 


f ind ] 
~ cos"6 
2t0 
cos @ | 


=Rceos™ 164 








cos” 6-cos né-isinn@ 





—isin@ 


sinné 


sin@ 


Also solved by D. A. Breault, Waltham, Massachusetts; Brother T. 
Brendan, St. Mary’s College, California; Huseyin Demir, Kandilli, Eregli, 
Kdz., Turkey; Joseph D. E. Konhauser, State College, Pennsylvania; 
Michael J. Pascual, Burbank, California; A. Kk. Rajagopal, Indian Institute 
of Science, Bangalore, India; Robert J. Wagner (Two solutions), Lebanon 
Valley College, Pennsylvania; Chih-yi Wang, University of Minnesota and 


the proposer. 


QUICKIES 


From time to time this department will publish problems which may be solved by 
laborious methods, but which with the proper insight may be disposed of with dispatch. 
Readers are urged to submit their favorite problems of this type, together with the elegant 


solution and the source, if known. 
gn 
Q 220. Show that lim S 1/k=log(q/p) (Submitted by Joseph Andrushkiw] 


N00 
Fa 
z a x 


] 
Q 221. Find L(x) whereit is defined by L(z) = lim -(sinzcos=+2sin=cos + 
Noo N n 2 ni 
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e 
don nsin— cos 2), [Submitted by Barney Bissinger] 


Q.222. On a 26-question test, 5 points were deducted for each wrong answer 
and § points were credited for each correct answer. If all questions were 
answered, how many wefe correct if the score was zero? [Submitted by 
C.W. Trigg] 


" esing 
Q223. Evaluate If 7 dx [Submitted by M.S. Klamkin} 
og 1+ cos*2 
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THE JAPANESE ABACUS 
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